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ABSTRACT 

The changes of hadron properties in dense and/or hot matter are discussed in 
terms of effective chiral Lagrangians with the parameters of the theory scaled in 
a simple way. The phenomenologically successful Walecka model is identified as 
a mean field chiral Lagrangian with the scaled parameters. Kaon condensation 
and chiral restoration transitions can be described within the same mean field 
framework. 
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Quantum chromodynamics (QCD) tells us that most, if not all, of light hadron masses 
are generated spontaneously by the breaking of chiral symmetry from SU(Nf) x SU(Nf) 
to diagonal SU(Nf) where Nf is the number of flavors, equal to 2 without strangeness and 
3 with. It is also widely believed that as a hadronic system is heated to high temperature 
or compressed to high density, the broken symmetry will get restored in a way paralleling 
what happens in condensed matter physics. A natural consequence of the restoration of 
the chiral symmetry must then be that the spontaneously generated masses disappear as 
density (and perhaps also temperature) is increased. The question we are raising is how 
does this "shedding of mass" occur? 

This question is at the core of the fundamental theory of matter: How is the mass 
generated, starting with the lightest object like neutrinos to the heaviest detected particle 
like the top quark? 

The aim of this talk is to describe how hadron properties get modified in medium as the 
system is heated or compressed. That is, immerse a hadron in medium and compress the 
system or heat it. What does one expect to see happening? 

To answer this question, let me start with the simplest nuclear system, namely the 
deuterium. Let us look at what happens when a soft photon is sent in to probe the system. 
Consider therefore the well-known inverse process 

n +p — > d + 7 (1) 

at thermal neutron energy. This process was first explained in a quantitative way by Riska 
and Brown [jij in 1973. What I will do here is to describe it more accurately in a modern 
QCD framework. Since it is a very low-energy process, QCD can be represented by an 
effective chiral Lagrangian field theory. This is because at long wavelength limit, chiral 
perturbation theory (ChPT) is believed to be exactly equivalent to QCD Q. This invites 
us to attempt to describe (|l]) in terms of a chiral Lagrangian. There is one basic problem 
in doing this and that has to do with the description of the deuteron in QCD: we do not 
really know how to derive the deuteron starting directly from QCD. Since we are focusing 
on the chiral aspect of the problem, however, the solution might be sought in a chiral 
Lagrangian approach to the structure of the deuteron. If the number of colors N c is in some 
sense big - which gives rise to what is known as the "large N c limit" - then the effective 
Lagrangian is given by meson fields only as we know from the skyrmion structure. There 
is some important progress in obtaining a bound deuteron, recently through the work of 
Manton |], as a baryon number 2 skyrmion, but we are still far from understanding it 
quantitatively. However there is an indirect approach to this which is consistent with QCD 
and which has the potential to be quantitatively accurate, namely that as we have known 
all along, the deuteron is made up of a proton and a neutron bound by meson exchanges: 
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In the framework of QCD, the nucleon may be gotten from a large N c Lagrangian as a 
soliton (skyrmion) but this is now known to be equivalent, at least in the large N c limit, to 
having the nucleon as a matter field in the chiral Lagrangian. We are thus led to consider 
a chiral Lagrangian that contains baryons (nucleon, A etc.), pseudo-Goldstone bosons tt 1 
(pions, kaons etc.), vectors = u^, p M , • • • etc. with suitable chiral invariant couplings. 
There have been some attempts to compute the deuteron from such a Lagrangian in chiral 
perturbation theory Q but at present the calculation can be done only at low orders since 
higher order calculations would involve too many parameters to be completely determined 
from available experiments. Luckily for our purpose, we need not compute the deuteron 
from first principles as I shall argue below. 

The important point to note is that certain aspects of the deuteron which have to do 
with the chiral symmetry structure of hadrons can be probed by the process (|l|) without 
knowing how to get the nucleus itself from a chiral perturbation theory, as recently discussed 
by Park, Min and Rho ||. Briefly the argument goes as follows: For physics with energy 
scale much less than the chiral scale A x ~ 4^/^ ~ 1 GeV (where f n is the pion decay 
constant ~ 93 MeV), the relevant Lagrangian is, schematically, 

C = Y,£i[B,U,M] (2) 

i 

where Bi are the baryon fields (both octet and decuplet), U the unitary pseudo-Goldstone 
fields U = e m l$T and A4 the quark mass matrix. The degrees of freedom more massive 
than the chiral scale A x are integrated out, appearing implicitly in the counter terms of the 
Lagrangian. Chiral symmetry requires that there be only derivative couplings apart from 
terms involving the mass matrix and hence effectively the Lagrangian is an expansion in 
d/A x and A^/A?,. Since the baryon mass is ~ A x , the baryons should be introduced as 
static matter fields, so the derivative on the baryon field does not involve time derivatives. 
ChPT is then a systematic perturbative expansion in powers of Q, say, Q n , where Q is the 
energy-momentum scale being probed, with suitable counter terms to remove divergences 
and to take account of the degrees of freedom that are integrated out. Now in applying this 
theory to nuclear systems, we need to separate the class of Feynman diagrams into two, one 
"irreducible" and the other "reducible." It is in calculating the irreducible diagrams that 
ChPT enters. The reducible diagrams - that cannot be treated by ChPT because of infrared 
singularity - are incorporated by solving a Lippman-Schwinger equation or Schrodinger 
equation with the potential obtained with the irreducible graphs by ChPT. This is how 
bound states are to be treated in ChPT. Now in calculating the process (jl]), we can write 
the EM current in two terms, one the single-particle current j' 1 ' and the other the two-body 
current 
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The former is called "impulse approximation current" and the latter "exchange current." 
For the system considered, we terminate with the two-body current. Later we will see that 
in heavy nuclei there can enter many-body currents, some of which become quite important. 

Very accurate wave functions for the final deuteron and the initial neutron-proton system 
obtained from some accurate phenomenological potential such as the Argonne v\s potential 
H (or a potential calculated in a high order ChPT if it is feasible) would correspond to 
a high order chiral expansion since the Schrodinger equation sums a certain class of chiral 
series to all orders and presumably the phenomenological potential also subsumes all orders 
of chiral perturbation. Now the idea is to compute the matrix element of the current (||) in 
chiral perturbation expansion in such a way that is consistent with the calculation of the 
wave functions. The calculation of the one-body current i« is without ambiguity. That of 
the two-body exchange current is somewhat subtle, requiring a careful sorting of irreducible 
and reducible contributions such that the reducible ones are suitably accounted for in the 
one-body term with the accurate wave functions. In this way, the Ward identities associated 
with the conserved vector current are satisfied to a given chiral order in the EM current. 
This is a numerically accurate procedure. 

A recent calculation || of the process (|l|) is given in Fig|l|. The calculation was done 
to order Q 3 which corresponds to next-to-next-to leading order in chiral expansion, that 
is, to one-loop order. Given the nucleon mass ijin, the pion decay constant f n , the pion 
mass m n , the axial coupling constant qa and the vector meson mass my all determined in 
free space, all the parameters that appear in the theory are fixed in the theory except for 
the hard core radius r c in the wave function reflecting on our inability to handle in ChPT 
very short-range physics. The remarkable agreement with the experiment shows that the 
two-nucleon systems we are looking at are made up of two nucleons with their properties 
as given in free space: the chiral Lagrangian with the vacuum values of hadron parameters 
describes nature remarkably well at low chiral orders. 

A very valuable information for heavier and denser nuclei is lodged in some of the terms 
that are negligible in the process ([]]) . As I showed elsewhere Q , two-body currents involving 
four baryon fields (e.g. Fig.||(c) with the blob replaced by a point) like 

BT^BBTB (4) 

with r and representing some Lorentz scalar and Lorentz vector quantities consistent 
with chiral symmetry, respectively, are subdominant in the chiral counting and can be 
ignored. Now terms like (Q) show up in the chiral Lagrangian as a result of integrating out 
heavy degrees of freedom with an energy scale E > A x . Consider for instance the scalar 
meson a that plays an important role in effective field theory of nuclei, a prototype of which 
being the Walecka model. In free space, there is no low-lying scalar that can appear in the 
low-energy chiral Lagrangian. But there is a high-lying scalar field that can be associated 
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Figure 1: The np capture rate calculated in chiral perturbation theory. The predicted 
capture rate (upper plot) a = (334 ±3) mb agrees with the experimental value (334.2 ±0.5) 
mb. The lower plot shows the ratios of the matrix elements of the twodoody terms over 
the one-body term. "Tree" corresponds to the leading chiral order one-pion exchange term 
(with the blob in Figs.||(a,b) replaced by a bare vertex), "l-7r(a;)" and "1tt(A)" correspond to 
the next-to-leading order corrections involving the uj meson and A resonance, respectively, 
to the one-pion exchange tree terms Figs.|2|(a,b). "27r" is the genuine loop correction to the 
tree contribution. The hadrons appearing in this calculation all have free-space properties. 
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Figure 2: Generic graphs contributing to exchange currents, (a) and (b) are one-pion 
exchange and (c) represents multipion and/or heavy- meson exchange currents. The large 
filled circles represent one-nucleon, one-pion irreducible graphs, the solid line the nucleon, 
the dotted line the pion and the wiggly line the current j M . 



with the trace anomaly of QCD 

(TpQCD 



/3(g) 
2ff 



fJ,V 



> 2 ~x 4 



(5) 



where (3 is the QCD beta function, g the color gauge coupling constant, G a ^ v the gluon field 
tensor and x t ne scalar glueball field. Now the \ field is massive, with m x ~ 2 GeV, so this 
degree of freedom appears only in the counter terms. It will give rise to a term like 



KBjjpBBB 



(6) 



with the coefficient k suppressed by the power (Q/m x ) 2 . In the process (||), Q is of order 
of 40 MeV, so (Q/m x ) 2 < 1CT 3 . Stated differently an effective two-body term like (||) will 
be screened by the short-range correlation implicit in the wave functions. 

The situation is quite different, however, in dense nuclear medium. As discussed in ||, 
as density increases, the scalar (or precisely the quarkish component of the scalar) x moves 
downwards in energy and at some high density, it joins the triplet of nearly massless pions 
to make up the quartet of the 0(4) symmetry of chiral symmetry. The merging presumably 
takes place at the chiral transition point discussed below. The point is that as discussed by 
Beane and van Kolck ||], in order to reconcile Weinberg's "mended symmetry" |l(]] with 
effective chiral Lagrangians at some shorter length scale, the scalar field must come down as 
a dilaton. Now if it comes down below the chiral scale Ay, then we can no longer consider 
the scalar as a counter-term contribution. It will strongly couple to low-mass multipion 
excitations giving among others what is usually taken as a scalar field a in effective nuclear 
forces ||. All other hadrons (other than Goldstone bosons) will also couple to this scalar as 
well and will undergo a mass shift as density increases. A chiral Lagrangian that accounts 
for this phenomenon has been shown to lead to the Brown-Rho scaling flT|| 



m* B /m B ~ my I my « m*/m a sa f*/f n 



(7) 
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where B stands for baryons, V for vector mesons and a for the dilatonic scalar. The star 
stands for density-dependent quantities. Such a chiral Lagrangian effective in dense system 
will then contain these effective constants instead of the free-space values used in the process 
© while preserving the free-space chiral symmetry. As I will argue, the consequence of this 
scaling can be significant in heavy nuclei and nuclear matter. 

One can already see the effect of this scaling in finite nuclei. One clear case is the axial 
charge transition in nuclei. It was shown in [l^] that to order Q 3 in chiral expansion the 
axial charge transition matrix element in heavy nuclei 

4(0+) ^ 5(0^), AT=1 (8) 

is enhanced with respect to the impulse approximation by the factor 

e M EC = ^(l + R) (9) 
m* N 

where R is the ratio of the exchange current matrix element Fig.^ to the impulse approxima- 
tion calculated with matter- free-space constants. In (|8]), R is essentially given by Figj2|(a) 
with the bare coupling and is given by R ~ 0.5 with a small variation with density of the 
system. Now at nuclear matter density, m^/uiN ~ 0.75, so 

e M EC « 2. (10) 
In light nuclei, we expect 6mec ^ 1-5. These results are in agreement what was found 



experimentally [13]. 



One predicts a similar effect in magnetic moments of heavy nuclei but here one has to 
include other effects of equal importance present in the vector current case. For instance, the 
"back-flow" correction due to Galilean invariance cancels almost completely the corrections 
coming from the scaled nucleoli mass. 

We now turn to an important issue of making a bridge between the chiral theory and 
Walecka mean field theory of nuclei and nuclear matter which is found to be very success- 



ful. I wish to show here, following [14], that Walecka theory is equivalent to the chiral 
Lagrangian theory at mean field with the BR scaling and that this would allow a treat- 
ment of fluctuations into different flavor directions (such as strangeness) in a way consistent 
with the properties of normal nuclear matter. To do this we can focus on the four-fermi 
interactions allowed in chiral Lagrangians that are relevant in making contact with Walecka 
theory, 

C Af = a(BB) 2 + f3(Bv»B) 2 (11) 

where a and (3 are dimension -2 constants and we are using the heavy-fermion formalism so 
that Vy, is the velocity four-vector of the heavy baryon. Let us imagine that the first term 
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of (JTT|) arises from integrating out the heavy chiral singlet scalar x anci the second term 
from integrating out a heavy chiral singlet vector meson uj. We can include other degrees 
of freedom in a similar way but we will not need them for symmetric nuclear matter that 
we shall consider. In this case, we can identify the constants 

Q 2 a 2 
« = A-' / 3 = -^TT- ( 12 ) 

We now ask what happens to this Lagrangian when it is immersed in dense and/or hot 
matter. In mean field, we get the nucleon scalar potential Sn and vector potential Vjv as 

9* 2 

Vn = (14) 



where p s is the scalar density and p the vector density. In obtaining (14), we have used 
SU(3)f relations together with KSRF relation which is known to hold well and put the 
stars in (|l<|) and (|i"4] ) to indicate that they are in-medium quantities. Now comparing with 
the phenomenology with Walecka model, we find that the identification requires that © 
holds with 

^-f r* 0.77. (15) 

Together with what we found in the case of the axial charge transition, we come to the 
conclusion that the mean field Walecka theory is just the mean field chiral Lagrangian 
theory with BR scaling. 

It is a well-known defect of the mean field Walecka model that the compressibility 
modulus Kq is much too high in the model. Now how does this defect get rectified? The 
answer must lie in higher loop corrections going beyond the mean field as the scaling is 
known to fail to give the nuclear matter saturation [Jl5||. This is also seen in recent work of 



Furnstahl et al [16] who note that by giving an anomalous dimension 2.7 to the scalar field a 
with the Lagrangian suitably implemented with the trace anomaly, they can obtain the low 
Kq ~ 200 MeV and the suppression of the many-body terms a n , n > 2. It is plausible that 
the anomalous dimension is mocking up the quantum loop effects that seem to be needed in 
the mean field approach given in 14]. A remarkable observation is that at the anomalous 
dimension of d a s» 2.7, two things happen simultaneously. One is that the Kq which is large 
at smaller anomalous dimensions stabilizes at ~ 200 MeV for d a ~ 2.7, stays at that value 
for higher <i a 's and secondly, it is at this fine-tuned value of d a that all multi-body forces 
get suppressed. This clearly calls for a simple explanation J!?]]. 

The above result immediately suggests how to calculate kaon-nuclear interactions in 
consistency with the nuclear matter properties as given by Walecka theory. To see this, 
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consider a part of the chiral Lagrangian that figures importantly in the kaon-nuclear sector 
Ckn = ^ (B lQ B)Kd t K + -jf (BB)KK = £ w + £ a (16) 

where K T = (K + K°). The constant / in (|l6|) can be identified in free space with the pion 
decay constant /„■. In medium, however, it can be modified as we shall see shortly. In chiral 
perturbation expansion, the first term corresponds to 0(Q) and the second term to 0(Q 2 ). 
There is one more 0(Q 2 ) term proportional to d 2 which will be taken into account in the 
numerical results quoted below but they are not important except for quantitative details. 

One can interpret the first term of ( |l6|) as arising from integrating out the u> meson as 
in the baryon sector. The resulting K~N vector potential in medium can then be deduced 
in the same way as for Vn- 

Vk± = ±^P- (17) 
Thus in medium, we may set / « /* and obtain 

V K ± = ±^V N . (18) 

This just says that the u couples to a matter field kaon, hence 1/3 of the ui coupling to 
the nucleon. The reason for this matter-field nature of the kaon is that all nonstrange 
hadrons become light in dense medium, so the kaon becomes in some sense heavy. This 
dual character is known from the hyperon structure which is well described by considering 
the kaon to be heavy as in the Callan-Klebanov model. 

As for the second term of (16), we use that the kaon behaves as a massive matter field. 
We therefore expect that it be coupled to the chiral scalar % as 

L a = hm K g*KK X (19) 

where the factor 1/3 accounts for one non-strange quark in the kaon as compared with three 
in the nucleon. When the \ field is integrated out as above, we will get, analogously to the 
nucleon case, 

1 q* 2 

C* = 2m K -^BBKK. (20) 
3 m* z 

Comparing with the second term of (|l6|), we find 

f 2 3 mf 1 ' 



We can get the T<kn from lattice calculations [18|, T,kn ~ S.2m n . This gives / w /*. 
Therefore we have 

S K± = 2 S N- (22) 
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To summarize: the kaon-nuclear potential gotten from a chiral Lagrangian and the 
nucleon-nuclear potential given by Walecka mean field theory are directly related through 
BR scaling. 

Given Walecka mean fields for nucleons, we can now calculate the corresponding mean- 
field potential for -fT~-nuclear interactions in symmetric nuclear matter. From the results 
obtained above, we have 

S K -+V K - ~~(S N -V N ). (23) 

Phenomenology in Walecka mean-field theory gives (SV — Vn) ^ —600 MeV for p = p . 
This leads to the prediction that at nuclear matter density 

S K - + V K - < -200 MeV. (24) 

This seems to be consistent with the result of the analysis in K-mesic atoms made by 
Friedman, Gal and Batty [19] who find attraction at p ~ 0.97po of 

S K - + Vk- = -200 ± 20 MeV. (25) 

An immediate consequence of this mean field description of the kaonic sector is that 
kaons will condense in dense neutron star (or nuclear star) matter at a density 

Pc ~ 2p (26) 

as found by Lee et al [20] in ChPT to one- loop order. 

To conclude, I make a few remarks on the nature of the scaling properties. 

1. There are two points to the issue. One is that some or all light hadrons may be 
undergoing a downward mass shift as density or temperature is increased. The second 
is that the scaling is like (0). These are basically two different issues. To the extent 
that hadron masses are generated spontaneously, it is inevitable that at least some 
masses should drop. Indeed QCD sum rule calculations do predict the drop for the 
vector- meson mass pit], the most recent value being |22|] 



m*/m p = 0.78 ± 0.08, p « p . (27) 

On the other hand, the "universal scaling" (]7|), obtained at mean field, may not be 
strictly valid. In fact, large N c arguments imply m^/mjv ~ \J g\l 'gA{ftl '/*■)• K we 
look at the strict N c = oo limit, g\ = gA as we showed in [11], the scaling (Q) holds 
but we know that in nuclei, g\ ~ 1 and thus in finite nuclei we expect that the nucleon 
scales somewhat faster than the pion decay constant, at least up to p ~ pq. Whether 
the "universal scaling" (|7|) is consistent with nature remains to be seen. 
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Figure 3: The Li-Ko-Brown explanation of the dilepton data of the CERES collaboration. 
The dotted line is the theoretical prediction without the scaling of the p and to mesons and 
the solid curve with the BR scaling. The three lowest mass points are essentially given by 
Dalitz pairs, so the relevant data points are the ones for higher invariant masses. Note that 
the peak shown at M ~ 800 MeV is predominantly given by the u decay outside of the 
medium. 
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2. An intriguing question is how far the scaling (^) can be pushed in density and/or 
temperature. Can one use the mean field argument all the way to the chiral phase 
transition? 



This is a highly relevant question since there are arguments [23] that the second order 
chiral phase transition relevant to QCD with 2 flavors is of mean field type as in 3D 
Gross-Neveu model. This suggests that the mean field chiral theory or Walecka theory 
could be used to discuss the QCD chiral phase transition [24]. If correct, this theory 
will be very useful for studying phase transitions in heavy ion collisions. 

3. There are a large number of experimental projects to measure the mass shift of hadrons 
in dense matter, particularly at GSI, CEBAF and CERN using dilepton-pair produc- 



tion. Recent data from the CERES collaboration [25] on e + e~ pairs in S on Au 
collisions at 200 GeV/nucleon can be understood in terms of the scaled mass of the 
vector mesons in medium, primarily due to density effect. In a recent paper, Li, Ko 
and Brown 26] have shown that the enhancement of the produced lepton pairs ob- 
served in the range of invariant mass 300 MeV < M < 550 MeV can be explained 
simply by the BR scaling in the p mass as the pairs are produced mainly through 
7r + 7r~ — > p* — > e + e~. The fit is given in FigFJ. The analysis is a complex one involv- 
ing the assumption of an expanding fireball in chemical equilibrium, but the economy 
of the explanation and the quantitative success make it quite compelling. 
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